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Abstract 

The recent discovery of a lateral organization in cell membranes due to small struc- 
tures called 'rafts' has motivated a lot of biological and physico-chemical studies. A 
new experiment on a model system has shown a spectacular budding process with 
the expulsion of one or two rafts when one introduces proteins on the membrane. 
In this paper, we give a physical interpretation of the budding of the raft phase. 
An approach based on the energy of the system including the presence of proteins 
is used to derive a shape equation and to study possible instabilities. This model 
shows two different situations which are strongly dependent on the nature of the 
proteins: a regime of easy budding when the proteins are strongly coupled to the 
membrane and a regime of difficult budding. 
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1 Introduction 



Classical and over-simplified models of the cell reduces the membrane to a 
bilayer of lipids in a fluid state which is a solvent for the proteins of the 
membrane (|l|). But the cell membrane is a much more complex and inhomo- 
geneous system. The inhomogeneities come from a phase separation between 
small structures called 'rafts' and the surrounding liquid phase. These rafts 
have been discovered a decade ago and remain an important issue of cell biol- 
ogy but also immunology, virology, etc (y). A lot of biological studies concern 
the rafts and examine their composition , their in- vivo size (j^, their role in 
signaling (0) or in lipid traffic (0) for example. The raft is roughly a mixture 
of cholesterol and sphingolipid but the exact nature of the sphingolipid and 
its concentration can vary between different rafts. In any case and whatever 
its composition, the raft has different physical or chemical properties than the 
rest of the membrane. In this paper, we focus on this specificity which is at 
the origin of an elastic instability that we want to explain. 



Experimentally, the raft in vivo cannot be easily studied and artificial systems 
like GUV (giant unilamellar vesicle) (0) appear more appropriate. GUV con- 
sist in a membrane of lipids with the possibility of a raft inclusion. On these 
artificial systems, a better control of the experimental parameters can be ob- 
tained and explored. For example, they have been used to study the coupled 
effects of both the membrane composition and the temperature on the nucle- 
ation of rafts (j^. Recently, a new experiment on GUV with rafts has shown a 
spectacular budding process (jl^ induced by injection of proteins called PLA2 
(phospholipase A2). Before injection, the GUV membrane is in a stable, nearly 
spherical state. But, more precisely, high-quality pictures of vesicles reveal two 
spherical caps, one for each phase: the raft and the fluid phase (jllh . These 
two caps have a radius of the same order of magnitude (about 5 micrometers, 
depending on the experimental conditions) and are separated by a discontin- 
uous interface. Few seconds after injection of PLA2 with a micro-pipette in 
the vicinity of the raft, one observes a rather strong destabilization of the 
initial shape: the raft tries to rise. The discontinuity of slope at the interface 
between the two caps becomes more and more pronounced. This lifting can 
be strong enough to expel the raft from the vesicle. When two or three rafts 
are present initially, successive expulsions can be observed. We present here 
a theoretical treatment showing that the driving force of the deformation is 
the absorption of proteins which locally deforms the membrane. We neglect 
chemical reactions since we focus here on the early stages of the instability: 
the time-scale of the instability is small compared to the characteristic time 
of chemical effects. We restrict ourselves to the simplest model relevant for 
the experiment we want to describe. It involves standard physical concepts of 
membrane mechanics. The initial shape of the system is given by a minimum 
of the energy of the whole system (that is the inhomogeneous vesicle including 
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proteins). A linear perturbation treatment allows to examine the existence of 
another solution which may lead to a new minimum of energy. This approach 
is sufficient to predict the experimental observation of destabilization and to 
derive a concentration threshold for an elastic instability of the vesicle. The 
calculation presented in this paper concerns only the first stages of the insta- 
bility. Intermediate stages require at least a dynamical nonlinear calculation 
including possible chemical effects of proteins. The final stage can be achieved 
by a non-linear calculation or, in case of fission, by the energy evaluation of 
two separated homogeneous spheres following the strategy described in ([l^. 



Models of vesicles have been widely described in previous papers ()l4 Il5 

They vary depending on the physical interactions involved taken 
into account. The backbone of all models is based on the minimization of the 
average curvature energy of the bilayer, with the introduction of a possible 
local membrane asymmetry ([lil : E3i. A large number of shapes have been 
predicted in the past by this model (jlGl ). They suitably describe experimental 
results such as the various shapes of red blood cells. Other physical effects 
can be introduced, such as the difference of area between the two layers of 
the membrane Q.suggestmg a differential compressive stress in the bilayer. 
These effects are visible under suitable experimental conditions (jioi ). Here, 
our scope is to study quantitatively the protein-membrane interaction using a 



generalization of the Leibler's model (|2lf ) to an inhomogeneous system. It turns 



out that this model, which describes the proteins as defects on the membrane, 
leads to a spatially inhomogeneous spontaneous curvature which is shown to 
be responsible for the destabilization of the whole system. Going back to the 
microscopic level, we derive a threshold for the protein concentration, which 
appears as a control parameter. Moreover, depending on the shape of the 
proteins, we are able to select two different regimes: a protein-stocking regime 
and a destabilization regime with possible raft-ejection. The idea of a non- 
homogeneous spontaneous curvature is not new since it has been used for 
mono-phasic vesicles to explain a possible thermal budding (jj^ . This does 
not concern the experiment described in ([lol ) since the temperature is not 
the relevant control parameter. Another scenario for the budding process of 
a raft has been proposed by ([lih: the increase of line tension by the proteins 
leads to an apparent slope discontinuity and to a neck. Again, this approach, 
which is different from ours, is not quantitatively related to the amount of 
proteins. It is why we suggest a different treatment as an interpretation of the 
raft ejection. 



This paper is organized as follows. Section 2 is devoted to a detailed description 
of the model defining precisely the elastic energy plus the energy of interac- 
tion combined to the constraints. Section 3 determines an obvious solution 
of the minimization of energy in terms of two joined spherical caps. A linear 
perturbation is performed which gives the threshold of proteins when a desta- 
bilization occurs. In section 4, the results are analyzed and discussed taking 
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into account known or estimated orders of magnitude of physical parameters. 



2 The model 



2.1 Membrane description 



The energetic model of the membrane is well established nowadays. It can 
incorporate many different interactions, constraints or restrictions. Here, we 
focus on a precise experiment and we think the model suitable for this exper- 
iment (0). Nevertheless, it can be modified easily for another experimental 
set-up. 

We consider a slightly stretched vesicle made of amphiphilic molecules diffi- 
cult to solubilize in water. The raft will be denoted by phase 1, it is usually 
considered as an ordered liquid. The remaining part is denoted by phase 2 and 
is considered as a disordered liquid. As the two phases are liquid, we describe 
them by two similar free energies, each of them having its own set of physical 
constants. Quantities which remain fixed in the experiment are constraints 
expressed via Lagrange multipliers in the free energy. So we define the energy 
of the bilayer in the phase (i): 



dS 



with H the mean curvature and K the Gaussian curvature. The square of H 



is the classical elastic energy ()14( ) when we get rid of the spontaneous curva- 
ture. Here, there is no physical reason to introduce a spontaneous curvature, 
sign of asymmetry between the two layers. The membrane contains enough 
cholesterol, which has a fast rate of fiip-fiop and which relaxes the constraints 
inside the bilayer. As for the Gaussian curvature, when a bi-phasic system 
without topological changes is concerned, it gives (Gauss- Bonnet theorem) a 
mathematical contribution only at the boundary (3)- means the surface 
tension: it is the combination between the stretching energy of the membrane 
and an entropic effect due to invisible fluctuations (jj^ . In addition to this 
energy of the bare membrane, we need to introduce the protein-membrane 
interactions. 



2.2 Protein-membrane interactions 



Both phases absorb the proteins, as soon as they are introduced, but probably 
with a different affinity. These proteins are not soluble in water, so we think 
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that they remain locahzed on the membrane and neglect possible exchange 
with the surrounding bath. As a consequence, the number of these molecules 
remains constant. Moreover, we assume that the proteins can not cross the 
interface (j2J). As suggested by S. Leibler (|2l[). the average curvature is coupled 
to the protein concentration, for two possible reasons. First, this can be due 
to the conical shape of the proteins which locally make a deformation of the 
membrane. Second, an osmotic pressure on the membrane results from the 
part of the protein in the water 0. Whatever the microscopic effect, the 
proteins force the membrane to tilt nearby and thus induce a local curvature. 
We define the energy due to proteins in the phase {i): 



/ a- 

AiH(f)+Xi(f)+ I yl 



dS 



(2) 



with (f) the concentration of proteins on the surface. The coupling constant 
is Aj. In Eq.(2), Aj is a Lagrange multiplier which allows to maintain the 
number of proteins constant in each phase. The model can be easily changed 
by considering Aj as the chemical potential of the proteins. In this case, the 
proteins are free to move everywhere on the membrane, to cross the interfaces 
or to go in the surrounding water. The last term in Eq.(2) is a Landau's 
expansion of the energy needed to absorb proteins on the surface nearby the 
equilibrium concentration (peq- The gradient indicates a cost in energy to pay 
for a spatially inhomogeneous concentration. Since the two phases are coupled 
together to make a unique membrane, let us describe now the interaction 
between them. 



2.3 Two phases in interaction 



The total energy of this inhomogeneous system is the sum of these two individ- 
ual energies for each phase plus at least two coupling terms. First, a more or 
less sharp interface exists between the raft and the phospholipidic part of the 
membrane. The interface is a line, the cost of energy of the transition being 
given by a line tension a equivalent to the surface tension in a vapor-liquid 
mixing. Second, the surface of the vesicle is lightly porous to the water but 
not to the ions or big molecules present in the solution. So, the membrane is 
a semi-permeable surface and an osmotic pressure appears. As a small varia- 
tion of the composition of the medium surrounding the vesicle induces a large 
variation of the size of the vesicle, the volume does not change when the pro- 
teins are injected, however the membrane will break down or transient pores 
will appear (jioh. which is not observed here. So, one needs to introduce a La- 
grange's parameter —P to express the constraint on the volume. Physically, P 
is the difference of osmotic pressure between the two sides of the membrane. 
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Fig. 1. Parameterization of an axisymmetric vesicle with two phases in cylindrical 
coordinates. 

Then, the free energy of the bilayer becomes: 



with C the boundary between the two phases. A variational approach is used 
to find the initial state and to study its stability to small perturbations. 



3 Static Solution and Stability analysis 

3. 1 Initial state 

First, we look for the simplest realistic solution with an homogeneous con- 
centration of proteins in each phase. It can be found by minimizing the en- 
ergy Ftot- This minimization gives the Euler-Lagrange equations (E-L equa- 
tions) plus the boundary conditions in an arbitrary set of coordinates. The 
surface has initially an axis of symmetry. Thus, the cylindrical coordinates 
seem to be the best choice. The parameterization of the surface is done by the 
arc-length s alone (see Fig.l). The energy becomes: 




i=l,2 



(3) 




(4) 
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with the new parameters: 



^'i = ^i + -^<Plq, and \[ = \i- aicpeqi- (5) 
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Minimization with respect to small perturbations of the spatial coordinates 
(r and ip) and of the protein concentration leads to the E-L equations: 



,,, sin('0) cos(-0) Pr , ib' , Aj ,, 7 , , , 

= — ' - TT- cos^ - — cos ^/^ + — 0' + — sin , 6a 

i = y V^" - ^ sin(^)' + - Pr sin(7/^) - A^^t/^' 

+f 0^ + + A>, (6b) 

0" = -|^(sin(^) + V^V) + ^0 - + ^, (6c) 

r' = cos('0). (6d) 

The boundary conditions deduced from Eq.(4) at the junction (defined by si) 
between the two phases and the continuity of the radius r give: 



K^ii^'isi - e)r(si) + (ki + kgi) sin(?/;(si - e)) - Ai0(si - e)r(si) (7a) 
-K2i)'{si + e)r(si) - {k2 + sin(V^(si + e)) + A20(si + e)r(si) = 0, 
7(si-e)-7(si + e)+a = 0, (7b) 
A0'(si - e)r(si) - /520'(si + e)r(si) = 0. (7c) 



The boundary conditions are deduced from the bounds in the variational pro- 
cess. The observation of a shape discontinuity at the boundary between the 



two phases (jllf ) suggests a solution which exhibits such a discontinuity of the 
slope at s = Si. The angle ip between the surface and the radius axis is chosen 
discontinuous at the interface (see Fig.2). This allows a tilt of the surface at 
the boundary C. The simplest solution, strongly suggested by the experiment, 
is two spherical caps of radius i?i and i?2, one for each phase with a constant 
concentration 0j (see Fig.2, for clarity, the deformation of the raft is stronger 
than the real one in the initial state). The minimization shows that one must 
satisfy the following "bulk" conditions for each phase in order to have this 
solution: 



X[Ri = 2Ai - aiRi(j)i so RiXi = 2Ai (8a) 
2E;i?2 _ pr3 ^ 2Ai(j)iRi - ai(t)lR^i = 0- (8b) 

The boundary conditions give two other relations: 
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Fig. 2. Parameterization of the vesicle near its initial state. The deformation is larger 
than the real one in the initial state. 

Ai0i = A2<^2, (9a) 

R2 005(^2) - Ri cos{e,) = (9b) 

FRi sm(t^i) 

where 61 and 62 are the polar angles in each phase at the boundary (see Fig.2). 
Note that, for each phase, the couple of equations given by Eq.(8) derive the 
Lagrange parameters like Aj (equivalent to a chemical potential) and Ej (the 
tension) which are quantities not easy to measure experimentally. On the 
contrary, Eq.(9) give geometrical informations. These informations with the 
other constraints such as the continuity of the radius and the ratio between 
area of both phases are enough to fix completely the values of -Ri, R2, Oi and 
^2- Both conditions have to be satisfied for all protein concentrations in order 
to ensure the existence of the initial homogeneous spherical caps, whether they 
are stable or not. If there is no protein, the conditions (8) reduce to A' = and 
2S-i?Q = P-Rq, which is the classical Laplace equation for an interface with a 
surface tension. Contrary to the law of capillarity, where the surface tension is 
a physical parameter dependent on the chemical phases involved, the tension 
here is not a constant characteristic of the lipids of the vesicle. It is a stress 
(times a length) which varies with the pressure. 



3.2 Linear perturbation analysis 



Now, we examine the stability of this solution. Due to the geometry, it turns 
out that the spherical coordinate system is more appropriate here and make 
the calculations easier. A perturbation of the spherical cap (i) is described 
by: R{6;(f) = Ri{l + u{6;(f)) with Ri the initial radius; in a similar way, a 
perturbation of the protein concentration is (p = 0,(1 + v{9;ip)) with 0, the 
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initial and homogeneous concentration of proteins on the surface. We assume 
that the hne tension is not modified by the addition of proteins, at least 
linearly. One can expand the free energy Eq. (4) to second order in u and v in 
the phase i: 



2tt 



/ Ci{u,'Vu, Au,v,'Vv) sinOdO 

J so, 



1 A 2 A 

-Au + -An^ + uAu H 

4 2 



-2A,;(/),;i?,; U 



Au Vu^ 



UV 



with 



(10) 



E^i?2 \2u + u' + 



Vu' 



f Am' 



PR" 



rv ■ 



^2/o.. , o.. , ..2 , ..... ,..2 , ^ ^^^2(^^)2 



2u + 2v + u^ + Auv + f ^ + 



Vu' 



3 
2 



(3n + 3^2 



+\'i(t)iR^i [2u + v + u^ + 2uv 



Vu' 



The total free energy is then a function of m and v, which allows a varia- 
tional approach to find the Euler-Lagrange's equations (E-L equations) and 
the boundary conditions. The E-L equations give shapes which are extrema of 
the free energy. Two sets of equations are derived: one for the zeroth order in u 
and V and one for the first order, the energy being calculated up to the second 
order of the perturbation. The zero-order equations gives the same results as 
Eq.(9). The first order equations are: 



Ai(l)iRi{2u + Au) + 0,2/?,' UiV - ^Av 







:ila) 



u 



2k, - E^i?2 + 2A,i?,0i - y 0,2/?,' - A>ii?2 ) An + k^AAm = 
equivalent to 



(lib) 



Ai(l)iRi{2v + Av) + [Ai^iRi 



PRl 



{2u + Au) + KiA {2u + Au) = (11c) 



This coupling imposes boundary conditions which must be treated at the zero- 
order and the first order. We have already studied the zero-order which gives 
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relations (9). As usual for linear perturbation analysis, the boundary condi- 
tions for the perturbation are homogeneous: u = Au = f = at the boundary 
between the two phases. Contrary to first intuition and usual procedures, al- 
though Eq.(lla) and (11c) are hnear, we cannot use the Legendre polynomial 
basis, due to the specific boundary conditions in this problem. The conve- 
nient angular basis in this case turns out to be the spherical cap harmonics, 
following standard techniques in geophysics (jj?! ) (see Appendix A where we 
recall some mathematical useful relations). These spherical cap harmonics are 
Legendre functions Pa;;(cos6'). The regular function at the pole of the cap is 
of the first kind and since we restrict on axisymmetric perturbations, these 
Legendre functions are simply hypergeometric function 



Notice that xi is not an integer. In the case where it is, we recover the Legendre 
polynomial basis. We select the spherical cap harmonics which vanish at the 
boundary angle [6 = 6i or 6 = 62, see Fig.2). This condition at the boundary 
gives a discrete infinite set of non-integer x\^^ values for the phase (i). I is an 
integer index used to order the allowed values x\^^ by increasing values. It is also 
the number of zero of the function P {») on the cap. These harmonics have the 



properties to be an orthogonal basis and to be eigenfunctions of the Legendre 
equation with eigenvalues: a:f^(xf'* + 1). We define u{6) = J2iUiPx,{cos6) and 
v{9) = T,iViPxi{cos 9) . From the first E-L equation (11a), one can deduce the 
amplitude Vi^i the protein concentration from ^ in the phase (i) : 



We introduce = x| (x; +l)/i?f , which is similar to the spatial period of the 
perturbation. Then, from the second E-L equation (lib) and after elimination 
using Eq.(12), we derive 



Our result can be compared to previous analysis made in two different asymp- 
totic limits in the homogeneous case. In these cases, the cap is a complete 
sphere and xi is an integer. First, for /3 = 0, we recover the result of ()28r i for 
an homogeneous vesicle without diffusion. Second, when Ri goes to infinity, 
we recover the result for an homogeneous fiat membrane (|21l). 

Notice that, in Eq.(13), the protein concentration has a similar significance 
as a negative surface tension: one can make the change of variable S^' = 




A, xr{xr + i)-2 






(13) 
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S' — aj0^/2 = Ej + aj(0egi ~ The principal effect of tlie proteins is to 

decrease tlie surface tension which is an obvious sign of instabihty. 



4 Discussion 

We will use the protein concentration (pi as our control parameter. Eq.(13) 
gives for each mode xf^ a threshold concentration $j such that for 0j < the 
initial state is stable and for 0i > one of the two phases is unstable, leading 
to a complete instability. The threshold concentration $j strongly depends on 
the physical properties of each phase. Then, the two parts of the initial system 
have no reason to be unstable simultaneously. The deformation of the other 
phase (not unstable to linear order) will be induced by the non-linear effects 
not included in this analysis. 

Since the thermal energy kT is the only external energy and the typical 
length of the phase (i) is -Rj, one can introduce dimensionless parameters: 



q = xjixi + l),^i = Ki/kT, t[ = J^iR^/kT, ~Ai = Ai/{kTRi), on = ai/{RjkT) 



and (3i = (3i/ [R^L^kT) with a characteristic length for the gradient of pro- 
tein concentration. Then, the protein concentration is replaced by Rfcff which 
is proportional to the square of the number of proteins in the phase (i). 

Rewriting the threshold (13), we find for the threshold concentration, in di- 
mensionless parameters: 



From Eq.(14), the search of the smallest threshold concentration gives two 
different regimes depending on the value of the dimensionless constant: 



This constant describes the strength of the coupling of the protein with the 
membrane (Aj) to the resistance of the membrane (kj) and to the absorption 
power (aj). In the weak interaction regime (c < 1), the protein concentra- 
tion is an increasing function of q (see fig. 3). So, the threshold is obtained 
for the smallest possible xf Xq. The direction of the deformation (inside or 
outside the initial cap) would be deduced from a third order calculation or 
from a numerical simulation. According to the definition of S- (Eq.(5)), the 
concentration required to destabilize the membrane is found bigger than the 
equilibrium concentration. But in this case, one expects that this threshold 




(14) 



c = 




(15) 
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Fig. 3. Threshold concentration versus the reduced modeg for two different values 
of (14) are taken as dj = 1, Ri = 1, t'- = 60, j3i{Lc/Rif = 0.01. The solid line is 
the case of large tilt due to each protein: Aj = 1.7. The dashed curve is the case of 
a smah tilt: Aj = 0.01. 

<l>j is difficult or impossible to reach since it requests the absorption of a con- 
centration of proteins larger than ^egr probably, the excess of proteins would 
prefer to dissolve in the surrounding water then forming aggregates. So, the 
weak regime of instability is not observable experimentally, our basic state 
made of two spherical caps is stable and proteins are stocked only. 

In the strong interaction regime, the limiting concentration shows a minimum 
not necessary for the smallest xi (see fig. 3). This has two consequences. First, 
the limiting concentration is less than the equilibrium concentration and it 
is easier to induce the instability. Second, the first unstable mode could be 
modified: g ~ 10 for the chosen numerical values. So, we have something 
more complex than the simple oblate/prolate (xq) deformation. This regime is 
observable and probably corresponds to the observed shape instability. In any 
case, our basic state cannot be seen in the experiment except as a transient. 

The existence of these two regimes, depending on the nature of the proteins, 
allows two possible and distinct scenarios for the cell: there is no doubt that 
this property is useful and probably used for biological purpose. The main 
difficulty of this study is the quantitative determination of the parameters 
since experimental values are not available even for this minimal model. Let 
us estimate c. The curvature of the membrane is of order 1/-Rj, its surface is 
close to R^. K is of order 10/cT for an unstretched vesicle, so k ^ 10. a, the 
cost of energy needed to increase the concentration of proteins, can be deduced 
from the energy required to remove a protein from the surface which is is about 
lOO/cT. It changes the concentration of proteins of 1/-Rf so 100/cT = ai/Rf 
and a ~ 100. (3 is deduced from which should be a small fraction of the 
radius of the sphere. We will take hereafter ~ -Rj/10. The proteins are 
moving at the surface of the membrane due to the Brownian motion. So the 
energy to move one protein by the length Lc is about kT. Then, Pi = kTL^Rf 
and /3i ~ 1. The value of Aj is more difficult to determine. Aj is the coupling 
constant between the membrane and the proteins. This is expressed by the 
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spontaneous curvature radius Rp. If Rp is small, the coupling effect is strong 
and, on the opposite, if Rp is large, Aj is small, which suggests that Aj is 
proportional to 1/Rp. But Aj is an energy multiplied by a length. Then, a 
good order of magnitude for Aj is kTRf/Rp, so Aj = Ri/Rp. Finally, we get 
c ~ Rp/Ri. li Rp is smaller than i?j, the system is in the strong coupling 
regime and in the other case, the interaction between membrane and proteins 
is weak. When the two phases have approximately the same physical constants 
( 111), the instability occurs first in the largest phase, as shown by Eq.(14) in the 
previous conditions. Nevertheless, the ejection of one part of the membrane 
requires a complete nonlinear dynamical treatment which will be derived from 
this energy formulation. 



5 Conclusions 



We have proposed a model of instability for an inhomogeneous vesicle which 
absorbs proteins. This instability is at the origin of a separation into two 
vesicles, one for each phase as seen experimentally. Our model rests on a 
"bulk" effect and assumes that the proteins are distributed everywhere on 
the membrane contrary to the "line tension "model which assumes a high 
concentration of the proteins at the raft boundary. To validate (or invalidate) 
our model, an experimental test could be the use of phosphorescent proteins 
with the same properties. It would be a way to follow the place where the 
proteins prefer to diffuse and stay on the membrane. Although we ignore the 
feasibility of such an experiment, it would provide a very useful information. 

We thank G. Staneva, M. Angelova and K. Koumanov for communicating 
their results prior to publication. We acknowledge enlightening discussions 
with J.B. Fournier. 



A Spherical cap harmonics 



The spherical cap harmonics are eigenvalues of the Laplace's equation in spher- 
ical coordinates. The Laplace's problem can then be rewritten as a Legendre's 
equation. The general solution is: 

t/r = /(0)L™(cos^) 

with 9 the colatitude, the longitude and an associated Legendre function. 
The eigenvalues associated to this solution are and —xi{xi + 1). So the 
solutions are symmetric with respect to xi = —1/2. So we can restrict to 
> —1/2 in all cases. 
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Generally speaking, m and xi can be integer, real or even complex and are 
determined by the boundary conditions. For a sphere, the solution must be 
periodic in the angle. This implies m real. In the particular case of an 
axisymmetric solution, which is the case in this paper, m = 0. 

The boundary condition on 6' for = is a condition of regularity: 



= for m = (A.l) 

oO 

t/^ = for m^O (A.2) 

It is satisfied by the Legendre functions of the first kind and excludes those 
of the second kind. Notice that this condition is required both for a complete 
sphere and for a spherical cap. 

In the case of the sphere, the boundary condition 6' = vr is similar to Eq.(A.l). 
The values of xi are then integer and the solutions are the classical associated 
Legendre polynomials. 

For a spherical cap whose ends are given by 6' = ±6'o, the boundary conditions 
at 6*0 are given by standard physical requirements. These boundary conditions 
can be satisfied by using two kinds of solutions such that either: 

^ = for d = ±do (A.3) 

or: 

^7^ = for 9 = ±00 (A.4) 
These conditions are satisfied by Legendre functions P^{cos9) with xi not 
necessary integer. No function can satisfy simultaneously the conditions (A.3) 
and (A.4) and there is two sets of xi which depend on the m value. We call 
yi{'m) the values of xi such as (A.3) is satisfied and zi{m) the values of xi such 
as (A.4) is satisfied. 

Functions in one set are orthogonal to each other but are not orthogonal to 
those of the other set. It is easy to show that: 



r P,7^(™)(cos^)P,7^(^)(cos^) sm9d9 = for h ^ k 

'J 

/ n7^(„)(cose)P™(^)(cose)sinerfe = for h^l, 

J 







P-(^)(cose)P-(^)(cos^)sin^cie = 

^oC, (»^) (cos ^) i [-P.; (cos 
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If the physics requires the boundary condition (A. 3) or (A. 4), the set of solu- 
tions Ui or Zi is enough to form a basis of solution of the problem. In the other 
case, one have to combine both of them and the resolution of the complete 
problem becomes more harder. In the case of this paper, we focus on the case 
of axisymmetric solutions (m = 0). The boundary conditions are given by A. 4, 
so the good set of parameters are the ?/i(0). The table A.l presents the first 
values of yi, calculated for two angles (tt/G and Stt/G), chosen as example. 
The figure A.l shows the three first P^, for 9q = tt/G. The figure A. 2 shows 
the three first P^^ for = Stt/G. The figure A. 3 shows the deformation of the 
spherical cap in the case of a perturbation by the three first P^;, for Oq = ir/G. 



I 





1 


2 


3 


4 


5 


6 


00 = vr/6 


4.08 


10.04 


16.03 


22.02 


28.01 


34.01 


40.01 


00 = 57r/6 


0.35 


1.57 


2.78 


3.98 


5.19 


6.39 


7.59 



Table A.l 

Values of yi for m = for the 7 first values of I. 



p 




Fig. A.l. Legendre functions Py^ versus the angle for = vr/6 and for / = (solid 
curve), I = 1 (large dashing) and / = 2 (small dashing) 



p 




Fig. A. 2. Legendre functions Py^ versus the angle for = 57r/6 and for I = 
(solid curve), I = 1 (large dashing) and I = 2 (small dashing) 
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Fig. A. 3. Deformations of a spherical cap of radius 1 and half-angle vr/G. The de- 
formations are dues to Legendre functions Py^ . For clarity, the maximum amplitude 
of the deformation is fixed to 0.1. This value is too large in principle for the linear 
analysis. The solid curve is the initial state. The largest dashing is the case / = 0. 
The intermediate dashing is the case / = 1. The smallest dashing is the case 1 = 2. 
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